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ABSTRACT. Recent phylogenetic studies of OpiHones have shown that Cyphophthalmi and Phalangida 
(= Palpatores + Laniatores) are sister groups, but higher relationships within Phalangida remain contro¬ 
versial. Current debate focuses on whether Palpatores (= Caddoidea + Phalangioidea + Ischyropsalidoidea 
+ Troguloidea) is monophyletic or paraphyletic, with Ischyropsalidoidea + Troguloidea (= Dyspnoi) being 
more closely related to Laniatores. The latter hypothesis was favored in recent combined studies of ri- 
bosomal DNA and morphology. Here higher relationships within Phalangida are examined using two 
nuclear protein-encoding genes, elongation factor-la (EF-la) and RNA polymerase II (Pol II), from 27 
opilion species representing seven super families. Cyphophthalmi was used as the outgroup. Nucleotide 
and inferred amino acid sequences were analyzed using maximum-parsimony and maximum-likelihood 
methods. All analyses recovered Palpatores as the monophyletic sister group to Laniatores with moderate 
to strong empirical support. Most palpatorean superfamilies were also recovered, but relationships among 
them were ambiguous or weakly supported. A monophyletic Palpatores was also obtained when EF~la 
and Pol II sequences were analyzed together with 18S and 28S rDNA sequences. 
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Members of the arachnid order Opiliones 
(harvestmen, shepherd spiders, daddy long- 
legs, etc.) are abundant and often highly vis¬ 
ible members of many terrestrial ecosystems. 
The group is estimated to encompass about 
5000 species (Shear 1982). The basic biology 
of most opilions is unexplored; but the order 
is known to encompass substantial behavioral 
and morphological diversity, with much recent 
work focusing on the structure and evolution 
of mating systems (e.g., Macias-Orddnez 
1997; Martens 1993; Mora 1987, 1990, 1991; 
Ramires & Giaretta 1994) and the evolution¬ 
ary morphology of genitalia (e.g., Martens 
1976, 1980, 1986; Martens, Hoheisel & Gotze 
1981; Hoheisel 1980; Shultz 1998). An un¬ 
derstanding of the phylogenetic relationships 
within Opiliones, and between Opiliones and 
other arachnid orders, is pivotal to further pro¬ 
gress in these and other fields of arachnolog- 
ical research. However, many relationships are 
unclear and the focus of ongoing investigation 
and debate (e.g., Giribet et al. 1995, 1999; Gi- 
ribet & Wheeler 1999; Shultz 1998, 2000). 


The goal of the present study is to reconstruct 
phylogenetic relationships among major line¬ 
ages within Opiliones using nucleotide and in¬ 
ferred amino-acid sequences from two nuclear 
protein-encoding genes, elongation factor-la 
(EF-la) and the largest subunit of RNA poly¬ 
merase II (Pol II). 

The monophyly of Opiliones is well estab¬ 
lished (Weygoldt & Paulus 1979; Shultz 
1990), and three principal opilion groups are 
widely recognized, namely, Cyphophthalmi, 
Palpatores and Laniatores (Hansen & Spren- 
sen 1904; Shear 1982). Cyphophthalmi (— Si- 
ronoidea) is a group of tiny (< 5 mm), hard- 
bodied and somewhat mite-like opilions 
characterized by many unique traits, such as 
a short male genital organ (spermatopositor), 
elevated ozopores (ozophores), a specialized 
spine on leg IV of males (adenostyle) and ab¬ 
sence of a genital operculum (Shear 1982; 
Shultz 1998). Laniatores is a species-rich 
group of heavily sclerotized opilions that have 
radiated extensively in the neotropics and 
southeastern Asia (Shear 1982). Members of 
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this undisputedly monophyletic group are 
characterized by large, often raptorial palps, 
unsegmented ovipositors and other unique 
features, Palpatores is a morphologically di¬ 
verse assemblage of four superfamilies (Cad- 
doidea, Ischyropsalidoidea, Phalangioidea, 
Troguloidea) united by few morphological 
synapomorphies„ Caddoidea and Phalangioi- 
dea include the typical “daddy longlegs” fa¬ 
miliar to inhabitants of northern temperate re¬ 
gions and are characterized, in part, by 
well-developed coxal lobes (coxapophyses) 
and segmented ovipositors. Ischyropsalidoi¬ 
dea and Troguloidea are less familiar groups 
but encompass a substantial range of morpho¬ 
logical diversity. They are characterized, in 
part, by reduced or absent coxapophyses, di¬ 
aphanous cheliceral teeth and unsegmented 
ovipositors. 

Opilion systematists have long debated the 
phylogenetic relationships of the three prin¬ 
cipal opilion lineages. Silhavy (1961) placed 
Cyphophthalmi and Palpatores together as the 
sister group to Laniatores on the basis of the 
arrangement of tarsal claws; that is, claws of 
all legs are similar in Cyphophthalmi and Pal¬ 
patores, but those of legs I and II differ from 
those of legs III and IV in Laniatores. A sim¬ 
ilar but more well-defended system was pro¬ 
posed by Martens and co-workers (Martens 
1976, 1980, 1986; Hoheisel 1980; Martens et 
ak 1981). They also united Cyphophthalmi 
and Palpatores (= Cyphopalpatores) but ar¬ 
gued that the palpatorean superfamily Trogu¬ 
loidea was the sister group to a clade com¬ 
prising Cyphophthalmi and the remaining 
Palpatores, a system that rendered Palpatores 
a paraphyletic group. Their Cyphopalpatores 
hypothesis was based largely on original mor¬ 
phological studies of the ovipositor and penis, 
character systems that are uniquely derived in 
Opiliones. Consequently, the root and branch¬ 
ing pattern of the tree proposed by Martens et 
al. was based on speculative character trans¬ 
formation series rather than inferences derived 
from outgroups. Shultz (1998) assessed the 
Cyphopalpatores concept through a parsimo¬ 
ny-based morphological analysis of genitalic 
and non-genitalic characters using a generic 
taxon sample similar to that of Martens et al. 
and outgroups (i.e., scorpions and xiphosu- 
rans) to polarize non-genitalic characters. The 
results supported Cyphophthalmi as the sister 
group to Laniatores + Palpatores (= Phalan¬ 


gida) Phalangida has also been supported by 
sperm ultrastructure (Juberthie & Manier 
1978), arrangement of extrinsic pharyngeal 
muscles (Shultz 2000) and, especially, molec¬ 
ular sequence data. Gi n bet et al. (1999) using 
18S and 28S ribosomal DNA and Shultz & 
Regier (unpublished) using EF-la and Pol II 
amino acid sequences have recovered both 
Cyphophthalmi and Phalangida as monophy¬ 
letic sister clades. 

Recent morphological and molecular anal¬ 
yses appear to have established Phalangida as 
a monophyletic group, but there is debate con¬ 
cerning basal relationships within this group. 
These opilions are traditionally divided into 
the Palpatores and Laniatores, with Palpatores 
often divided into Dyspnoi (= Ischyropsali¬ 
doidea + Troguloidea) and Eupnoi (= Cad¬ 
doidea + Phalangioidea) (Hansen & S0rensen 
1904; Silhavy 1961; Juberthie & Manier 
1978; Shultz 1998). However, Giribet et ah 
(1999) conducted a phylogenetic analysis of 
Opiliones using 18S and 28S ribosomal DNA 
sequences and recovered two topologies, one 
favoring a monophyletic Palpatores and one 
favoring a paraphyletic Palpatores, with Dy¬ 
spnoi being the sister group to Laniatores. The 
molecular characters alone did not strongly fa¬ 
vor one hypothesis over the other, but a mor¬ 
phology-based analysis and the combined mo¬ 
lecular-morphological study tended to support 
the Dyspnoi 4- Laniatores hypothesis. In a 
subsequent study, Giribet & Wheeler (1999) 
explored the significance of indels on the phy¬ 
logenetic analysis of Opiliones. Using 18S ri¬ 
bosomal DNA sequences and morphological 
characters from Giribet et al. (1999), Giribet 
& Wheeler recovered Dyspnoi + Laniatores. 
However, given their exclusion of the 28S 
rDNA and use of a problematic morphology 
matrix (see Methods for details), it is probably 
best to regard Giribet & Wheeler's contribu¬ 
tion as a demonstration of a new analytical 
method rather than a study of opilion relation¬ 
ships. 

The present study examines higher relation¬ 
ships within Phalangida using nucleotide and 
amino-acid sequences from EF-la and Pol II 
with the specific aim of testing the Dyspnoi 
+ Laniatores hypothesis. Sequences were ob¬ 
tained from 27 opilion taxa representing the 
major opilion superfamilies (i.e., Sironoidea, 
Travunioidea, Gonyleptoidea, Phalangioidea, 
Caddoidea, Ischyropsalidoidea and Troguloi- 
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dea). Maximum-parsimony and maximum- 
likelihood analyses strongly and consistently 
recovered Falpatores and Laniatores as mono- 
phyletic sister groups under all character par¬ 
titions, weighting schemes and analytical 
methods. This strong support for the mono- 
phyly of Falpatores contrasts with the ambig¬ 
uous or weakly supported conclusions derived 
from studies of 18S and 28S ribosomal DNA 
and morphology (Giribet et ah 1999; Giribet 
& Wheeler 1999) or morphology alone 
(Shultz 1998). The represented superfamilies 
except Travunioidea also tended to be recov¬ 
ered, although support for Ischyropsalidoidea 
was weak, and it was not possible to conclude 
whether Dyspnoi and Eupnoi are monophy- 
letic sister groups. Combined analysis of EF- 
la and Pol II sequences with the 18S and 28S 
rDNA sequences of Giribet et ah (1999) also 
consistently recovered a monophyletic Palpa- 
tores. We regard these results as compelling 
evidence in support of the Falpatores hypoth¬ 
esis and against the Dyspnoi + Laniatores hy¬ 
pothesis. We suggest that future work on high¬ 
er-level relationships within Opiliones focus 
an assessing the monophyly and relationships 
of superfamilies within Falpatores and Lan¬ 
iatores. 

METHODS 

Abbreviations.—Abbreviations used in the 
present study are as follows: bp, base pairs; 
ntl, first codon position; ntlnoLR, ntl data 
subset in which any characters that encode a 
leucine or arginine residue for any taxon are 
excluded at the homologous position in all 
taxa; ntlLR, ntl data subset which includes 
any characters that encode a leucine or argi¬ 
nine residue for any taxon plus all other ho¬ 
mologous characters in other taxa; nt2, second 
codon position; nt3, third codon position. 

Terminal taxa and sequences.—Sequenc¬ 
es of elongation factor-la (EF-la: 1092 bp) 
and the largest subunit of RNA polymerase II 
(Pol II: 1038 bp) were generated from 27 op- 
ilion species representing seven superfamilies. 
Specimens were collected alive, killed by im¬ 
mersion in 100% ethanol, and stored in 100% 
ethanol at —20 °C or — 80 c C. Voucher spec¬ 
imens will be deposited in the U.S. National 
Museum of Natural History (Smithsonian In¬ 
stitution, Washington, DC). The list of speci¬ 
mens follows with GenRank accession num¬ 
bers for EF-la and Pol II, respectively, in 


parentheses: Sironoidea: Parasiro coiffaiti 
Juberthie 1956 (Sironidae) (AF240852; 
AF241025 - AF241027), Siro acaroides (Ew¬ 
ing 1923) (Sironidae) (AF240855; AF241034 
- AF241036). Travunioidea: Equitius doriae 
Simon 1880 (Triaenonychidae) (AF240867; 
AF241068 - AF241070), Sclerobunus robus- 
tus (Packard 1877) (Triaenonychidae) 
(AF240858; AF241042 - AF241044). Gony- 
leptoidea: Gonyleptes fragilis Mello-Leitao 
1922 (Gonyleptidae: Gonyleptinae) 

(AF240868; AF241071 - AF241073), Pro- 
gonyleptoidellus striatus (Roewer 1913) 
(Gonyleptidae: Progonyleptoidellinae) 

(AF240873; AF241086 - AF241088), So- 
dreana sodreana (Mello-Leitao 1922) (Gon¬ 
yleptidae, Sodreaninae) (AF240879; 
AF241104 - AF241106), Promitobates orna- 
tus (Mello-Leitao 1922) (Gonyleptidae, Mi- 
tobatinae) (AF240876; AF241095 - 

AF241097), Discocyrtus areolatus Piza 1938 
(Gonyleptidae, Pachylinae) (AF240842; 
AF240994 - AF240996), Laneosoares inermis 
(R. Soares 1944) (Tricommatidae) 
(AF240871; AF241080 - AF241082), Pseu- 
dobiantes japonicus Hirst 1911 (Phalangodi- 
dae) (AF240874; AF241089 - AF241091), 
Proscotolemon sauteri Roewer 1916 (Phal- 
angodidae) (AF240872; AF241083 - 

AF241085), Scotolemon lespesi (Lucas 1860) 
(Phalangodidae) (AF24G878; AF241101 - 
AF241103). Caddoidea: Caddo agilis Ranks 
1892 (Caddidae) (AF240838; AF240980 - 
AF240982), Caddo pepperella Shear 1975 
(Caddidae) (AF240863; AF241056 - 

AF241058). Phalangioidea: Astrobunus gral¬ 
lot or Simon 1879 (Sclerosomatidae) 
(AF240862; AF241053 - AF241055), Odiel- 
lus picius (Wood 1868) (Phalangiidae) 
(AF240850; AF241019 - AF241021), Phal- 
angium opilio Linneaus 1758 (Phalangiidae) 
(AF240875; AF241092 - AF241094). Ischy¬ 
ropsalidoidea: Ceratolasma tricantha Good¬ 
night & Goodnight 1942 (Ceratolasmatidae) 
(AF240864; AF241059 - AF241061), Hespe- 
ronemastoma modestum (Ranks, 1894) (Cer¬ 
atolasmatidae) (AF240869; AF241074 - 
AF241076), Ischyropsalis luteipes Simon 
1872 (Ischyropsalididae) (AF240870; 
AF241077 - AF241079), Sabacon imamurai 
Suzuki 1964 (Sabaconidae) (AF240877; 
AF241098 - AF241100), Taracus pallipes 
Banks 1894 (Sabaconidae) (AF240881; 
AF241110 - AF241112). Troguloidea: Di- 
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cranolasma scabrum (Herbst 1799) (Dicran- 
olasmatidae) (AF240866; AF241065 - 

AF241067), Trogulus nepaeformis (Scopoli 
1763) (Trogulidae) (AF240880; AF241107 - 
AF241109), Nipponopsalis abei (Sato & Su¬ 
zuki 1939) (Nipponopsalididae) (AF137391; 
AF138993 - AF138995), Dendrolasma denti- 
palpe Shear & Gruber 1983 (Nemastomati- 
dae) (AF240865; AF241062 - AF241064). 

Outgroups and phylogenetic frame¬ 
work.—The present study focuses on higher- 
level relationships within Phalangida and thus 
uses Cyphophthalmi as the outgroup. As sum¬ 
marized in the introduction, Opiliones is an 
unambiguously monophyletic group and re¬ 
cent phylogenetic analyses strongly support 
Cyphophthalmi and Phalangida as monophy¬ 
letic sister groups. Results from molecular 
studies have been particularly convincing. Gi- 
ribet et al. (1999) conducted an analysis of 
higher relationships within Opiliones and out¬ 
groups (i.e., Ricinulei, Solifugae, Scorpiones, 
Xiphosura) using 18S and partial 28S ribo- 
somal DNA sequences and consistently recov¬ 
ered Cyphophthalmi and Phalangida as sister 
groups. Shultz & Regier (unpublished) used 
EF-la and Pol II in a broad study of ordinal 
relationships within Arachnida and recon¬ 
structed Cyphophthalmi and Phalangida as 
monophyletic clades with high bootstrap sup¬ 
port. Opiliones was rooted in the analysis with 
26 outgroup species representing Xiphosura 
and all arachnid orders except Palpigradi. Tak¬ 
en together, these results effectively falsify the 
Cyphopalpatores hypothesis and strongly sup¬ 
port the monophyly of Phalangida. It would 
have been possible in the present study to in¬ 
clude all or some non-opilion arachnids, but 
these were excluded because inclusion of cer¬ 
tain rapidly evolving arachnid lineages (es¬ 
pecially Acari and Pseudoscorpiones) desta¬ 
bilizes even well-established relationships 
within Opiliones. We have therefore used the 
two species from Sironoidea (Cyphophthalmi) 
to root Phalangida. 

Sequence data.—Procedures to generate 
sequence data sets have been published (Re¬ 
gier & Shultz 1997). In brief, total nucleic ac¬ 
ids were isolated; cDNA copies of EF-la and 
Pol II mRNA were reverse transcribed; ds- 
DNA copies were amplified by PCR and sub¬ 
sequently gel isolated; the resulting PCR frag¬ 
ments were used as templates for another 
round of PCR amplification with nested prim¬ 


ers; the resulting fragments were gel isolated 
and sequenced. If the resulting fragment con¬ 
centration was too low to sequence directly, it 
was either concentrated or reamplified with 
M13 sequences present at the 5' ends of all 
primer sequences. The same M13 sequences 
were also used as primers for thermal cycle/ 
dideoxy sequencing. Sequencing reactions 
were fractionated and preliminary analyses 
were performed with Perkin-Elmer/ABI au¬ 
tomated DNA sequencers. Automated DNA 
sequencer chromatograms were edited and 
contigs were assembled using the pregap and 
gap4 programs within the Staden software 
package (Staden et al. 1999). Sequences from 
different species were aligned and Nexus-for¬ 
matted nucleotide data sets were constructed 
using the Genetic Data Environment software 
package (version 2.2, Smith et al. 1994). All 
sequences lacked indels. Amino acid data sets 
were inferred from nucleotide sequences using 
the universal nuclear genetic code option in 
MacClade, version 3.08 (Maddison & Mad- 
dison 1992). 

Data analysis.—Maximum parsimony 
(MP) analyses of nucleotide and amino acid 
data sets were performed in PAUP*4.0 (Swof- 
ford 1998) using unordered, equally weighted 
characters. The following data sets were ana¬ 
lyzed: all-nucleotide, ntl + nt2, ntlnoLR + 
nt2, nt3, amino acids. Amino acids were also 
analyzed using a “Protpars” step matrix con¬ 
structed in MacClade. The step matrix con¬ 
sisted of transformation scores of “1” or “2” 
determined by the minimum number of non- 
synonymous nucleotide changes separating 
particular codons. Serine codons that differed 
at ntl were coded separately (termed “non- 
disjunct” in MacClade). Analysis consisted of 
a heuristic search using TBR branch swapping 
with random taxon addition (100 replications). 
Bootstrap values (Felsenstein 1985) were also 
obtained using heuristic searches (1000 rep¬ 
lications each with 10 random-sequence ad¬ 
dition replicates). The incongruence length 
difference test (Farris et al. 1995), imple¬ 
mented in PAUP* 4.0 as the partition homo¬ 
geneity test (1000 replications), was used to 
test for the significance of conflict between 
genes and character partitions. 

Maximum likelihood (ML) analyses of nu¬ 
cleotide data sets were performed with PAUP* 
4.0 under the optimal GTR model of sequence 
evolution, the optimal model based on a series 
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of likelihood-ratio tests conducted according 
to Swofford et al. (1996: 430-438; see also 
Shultz & Regier 2000). Among-site rate var¬ 
iation was accommodated within the GTR 
model by likelihood estimations of separate 
rates for individual codon positions by gene 
and by fitting total likelihood-estimated char¬ 
acter change to a gamma distribution with in¬ 
variable sites estimated separately (Hasegawa, 
Kishino & Yano 1985). We call the former 
model GTR-ssr and the latter GTR + T + I. 
The GTR-ssr model was applied to total nu¬ 
cleotide data and the GTR + F + I model to 
the ntl + nt2 and the ntlnoLR + nt2 data 
sets. 

As a first step in our exploration of tree 
space using ML analysis, we used an MP tree 
derived from analysis of amino acids as the 
input topology on which likelihood parame¬ 
ters were optimized. NNI branch swapping 
was then performed, and new likelihood pa¬ 
rameters were estimated from the most likely 
topology. TRR branch swapping was con¬ 
ducted on the new tree and likelihood param¬ 
eters were re-estimated. These parameters 
were then used as input for a heuristic search 
with NNI branch swapping and 100 random 
taxon additions. The parameters from the 
overall best tree were re-optimized. Bootstrap 
analyses were too computationally demanding 
to be performed in the same manner. Instead, 
the Neighbor Joining algorithm coupled with 
an ML-estimated distance matrix were used 
(1000 replications). For each bootstrapped 
data set, a ML distance matrix was calculated 
that assumed a minimum evolution objective 
function and that used identical parameters 
(i.e., rate matrix, base frequencies, a value, 
proportion of invariant sites) to those estimat¬ 
ed from the ML topology of the original data 
set. Bootstrap values for the ntl + nt2 and the 
ntlnoLR + nt2 data sets were calculated in 
this manner. 

ML analysis of the amino-acid data set was 
performed using the protml program within 
the MOLPHY software package (version 2.2, 
Adachi & Hasegawa 1994) and the empirical 
transition matrix compiled by Jones, Taylor & 
Thorton (1992). All 37681 amino acid parsi¬ 
mony trees within 4 steps of the minimum- 
length tree (= 354 steps) were read in batch 
into protml, and the tree with the highest like¬ 
lihood score was selected. The significance of 
differences in the fit of individual data sets to 


alternative topologies under both MP and ML 
criteria was assessed by the test of Kishino & 
Hasegawa (1989), using the “Tree Scores” 
option as implemented in PAUP* 4.0. Only 
fully dichotomous trees were compared; that 
is, single nodes (e.g., Dyspnoi + Laniatores) 
were constrained and the remaining set of un¬ 
constrained relationships was reoptimized. 
Percentage differences of all pairwise combi¬ 
nations of EF-la and of Pol II amino acid and 
nt3 data sets were calculated in PAUP* 4.0. 
Average differences across the basal node of 
various groups were calculated by averaging 
all values across the basal dichotomy within 
a particular clade. Base compositions were 
calculated by gene and by codon position type 
using PAUP* 4.0. 

Combined analyses with ribosomal 
DNA.—Unweighted MP analyses were con¬ 
ducted on data matrices that combined the 18S 
and 28S ribosomal sequences generated by 
Giribet et al. (1999) (15 opilions and 5 out¬ 
groups) with the EF-la and Pol II sequences 
(27 opilions) generated in the present study. 
Two combined matrices were constructed. In 
the first, ribosomal data were combined with 
complete nucleotide sequences from EF-la 
and Pol II for a total of 4329 characters (1117 
parsimony-informative characters), and, in the 
second, the ribosomal data were combined 
with amino acids from EF-la and Pol II for a 
total of 2909 characters (346 parsimony-in¬ 
formative characters). In both cases, the ri¬ 
bosomal sequences were aligned and edited as 
described in Giribet et al. (1999). Both matri¬ 
ces included 35 terminal taxa, which encom¬ 
passed 15 taxa with both ribosomal and EF- 
la + Pol II sequences, nine taxa with only 
ribosomal sequences, and 11 taxa with only 
EF-la + Pol II sequences. Regions of the ma¬ 
trix lacking sequence data were treated as 
missing characters. 

The specific strategy for combining ribo¬ 
somal and protein-encoding sequences is pre¬ 
sented below. The following five species were 
represented by both rDNA and protein-encod¬ 
ing sequences: Parasiro coiffaiti (GenBank 
accession nos.: 18S: U36999; 28S: U91495), 
Equitius doriae (18S: U37003, 28S: U91503), 
Scotolemon lespesi (18S: U37005, 28S: 
U91506), Caddo agilis (18S: U91487, 28S: 
U91502) and Ischyropsalis luteipes (18S: 
U37QG0, 28S: U91497). We also combined se¬ 
quences from five pairs of closely related taxa. 
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specifically, Siro rubens Latreille 1804 (18S: 
U36998, 28S: U91494) was combined with 
EF-la and Pol II sequences for Siro acaro- 
ides , Odiellus troguloides (Lucas 1847) (18S: 
X81441, 28S: U91500) was combined with 
Odiellus pictus , Dicranolasma soerenseni 
Thorell 1876 (18S: U37001, 28S: U91498) 
was combined with Dicranolasma scabrum , 
the nemastomatid Centetostoma dubium (Mel- 
lo-Leitao 1936) (18S: U37002, 28S: U91499) 
was combined with Dendrolasma dentipalpe , 
and the pachyline gonyleptid Pachyloides tho- 
relli Holmberg 1878 (18S: U37007, 28S: 
U91508) was combined with Discocyrtus ar- 
eolatus. Taxa represented only by rDNA se¬ 
quences included an unidentified Stylocellus 
Westwood 1874 (18S: U91485, 28S: 

U91496), an unidentified Oncopus Thorell 
1876 (18S: U91488, 28S: U91504), Maiore- 
rus randoi Rambla 1991 (18S: U37004, 28S: 
U91505), Gnidia holnbergi Soerensen 1912 
(U37006, 28S: U91507), the solifuge Eusi- 
monia wunderlichi Pieper 1977 (18S: 

U29492, 28S: none), the ricinuleid Pseudo- 
cellus pearsei (Chamberlin & Ivie 1938) (18S: 
U91489, 28S: none), the scorpion Androcton- 
us australis C.L. Koch 1839 (18S: X77908, 
28S: none), and the xiphosurans Limulus po- 
lyphemus (Linneaus 1758) (18S: U91490, 
28S: U91492) and Carcinoscorpius rotundi- 
cauda (Latrieille 1802) (18S: U91491, 28S: 
U91493). 

Giribet et al. (1999) had also presented a 
matrix of 45 morphological characters from 
an unreferenced literature review, but we did 
not include it in the combined analysis due to 
numerous ambiguities and inaccuracies. For 
example, palpal claws were coded as well de¬ 
veloped in Sironoidea, but these claws are re¬ 
duced compared to pedal claws (Hansen & 
Sprensen 1904; Shultz 1998). Two characters 
focusing on “fusion of abdominal tergites” 
were questionable, because males and females 
were coded separately, thereby making the 
characters non-independent. Metapeltidial 
cones were coded as “absent” in Caddo but 
are present (Shultz 1998). Internal longitudi¬ 
nal muscles of the ovipositor were coded as 
“present” in Parasiro and Siro ; “reduced” in 
Odiellus , Caddo , Dicranolasma , Pachyloides , 
Scotolemon and Centetostoma but they are 
“absent” in sironoids, phalangioids, Caddo 
and Dicranolasma and well developed in Cen¬ 
tetostoma , Scotolemon and other gonyleptoids 


(Martens et al. 1981). Circular muscles of the 
ovipositor were coded as absent in Odiellus 
and Caddo , but they are present in phalan¬ 
gioids and Caddo (Martens et al. 1981). Given 
concerns about accuracy, the absence of 
source citations and the need to expand the 
matrix to include taxa not considered by Gi¬ 
ribet et al. (1999), we chose not to include the 
morphological matrix in the combined study. 

RESULTS 

Pairwise differences and base composi¬ 
tion.—Inferences of phylogenetic relation¬ 
ships among arthropods at taxonomic levels 
deeper than those in the present study have 
demonstrated that EF-1 a and Pol II sequences 
retain phylogenetic signal, although the syn¬ 
onymous changes of nt3 can be very homo- 
plasious due to overlapping substitutions. 
Across Opiliones, observed pairwise differ¬ 
ences at nt2, at which all changes are non- 
synonymous, did not exceed 8.1% for either 
gene. By contrast, maximum observed pair¬ 
wise differences at nt3 were 50.6% for EF-la 
and 70.5% for Pol II, and were never less than 
10% for either gene. These observations are 
consistent with the greater homoplasy at nt3 
relative to nt2 (and nt 1, in which changes are 
a mixture of synonymous and non-synony- 
mous changes). For example, fitting ntl, nt2, 
and nt3 data sets to the topology shown in Fig. 
1 resulted in overall Retention Indices of 
0.6347, 0.8351, and 0.3825, respectively. 

Chi-square tests of homogeneity across taxa 
revealed no significant heterogeneity (P » 
0.05) in most partitions examined (i.e., ntl, 
nt2 and nt lnoLR of both genes separately and 
combined), both with and without constant 
sites excluded. Exceptions were ntl minus 
constant sites for Pol II, nt3 for Pol II, and 
those partitions that included nt3 of Pol II 
(i.e., all-nucleotides and nt3 combined), all of 
which showed highly significant base hetero¬ 
geneity (P < 0.0001). 

Phylogenetic analyses of EF-la and Pol 
II. —Partition homogeneity tests of character 
partitions revealed no significant inconsisten¬ 
cies in the phylogenetic signals of EF-la and 
Pol II, except in comparisons that included nt3 
of Pol II (Table 1). Given these results, and 
uncertainties concerning the importance and 
utility of the partition homogeneity tests as an 
assay of combinability (Baker & DeSalle 
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1995; Cannatella et al. 1998), we performed 
combined analyses of all partitions. 

Maximum-parsimony (MP) analyses of all 
character partitions (all-nucleotides, ntl +2, 
ntlnoLR + nt2, nt3, amino acids, amino-ac- 
ids-protpars) recovered the following higher 
clades within their respective minimal-length 
trees (bootstrap percentages [BP] in parenthe¬ 
ses): Phalangida (100, 100, 100, 95, 100, 
100); Laniatores (92, 99, 100, 45, 100, 100); 
Palpatores (100, 96, 79, 88, 87, 76); Trogu- 
loidea (93, 97, 94, 45, 94, 97); Caddoidea 
(100, 100, 100, 100, 100, 100) and Phalan- 
gioidea (100, 100, 100, 98, 100, 100). Ischy- 
ropsalidoidea was recovered by all-nucleo- 
tides (BP 56), ntl + nt2 (BP 73), ntlnoLR + 
nt2 (BP 43) and amino-acids-protpars (BP 30) 
but was recovered as a paraphyletic group by 
nt3 and amino acids. Dyspnoi (= Troguloidea 
4- Ischyropsalidoidea) was recovered by all¬ 
nucleotides (BP 24), ntl + nt2, (BP 59), 
ntlnoLR + nt2 (BP 54), and amino-acid-prot- 
pars (BP 47). Eupnoi (= Caddoidea + Phal- 
angioidea) was recovered by ntlnoLR +nt2 
(BP 36) within a subset of 20 MP trees. Gon- 
yleptoidea was recovered by all-nucleotides, 
ntlnoLR + nt2, amino acids, and amino-ac¬ 
ids-protpars with weak-to-moderate bootstrap 
support (i.e., 29-74%). Travunioidea was re¬ 
covered as a paraphyletic group with Equitius 
being the sister group to Gonyleptoidea and/ 
or occurring within Gonyleptoidea. 

The all-nucleotide and ntl + nt2 data sets 
significantly preferred a monophyletic Palpa¬ 
tores over the MP clade constrained to a 
monophyletic Laniatores + Dyspnoi (P values 
ranged from 0.015 to 0.028), according to the 
test of Kishino & Hasegawa (1989); the other 
data sets were indecisive. 

The ML topology recovered by analysis of 
the ntlnoLR + nt2 data set is shown in Fig. 
1. All maximum likelihood analyses (three 
models of nucleotide substitution and four 
data subsets, see Methods) recovered the fol¬ 
lowing clades in their ML topologies (BP val¬ 
ues for analysis of ntl + nt2 and ntlnlLR + 
nt2, respectively, are in parentheses): Phalan¬ 
gida (100, 100), Laniatores (100, 100), Pal¬ 
patores (96, 78), Troguloidea (98, 93), Gon¬ 
yleptoidea (63, 85), Phalangioidea (100, 100), 
Caddoidea (100, 100), and Ischyropsalidoidea 
(53, 42). Eupnoi (37, 55) was recovered by all 
three nucleotide data sets. Dyspnoi (19, 36) 
was recovered only by the ntlnoLR + nt2 


data set. ML analysis of amino acids recov¬ 
ered Phalangida, Laniatores, Palpatores, Tro¬ 
guloidea, Caddoidea, Phalangioidea, and Gon¬ 
yleptoidea. 

Combined analyses.—Combined un¬ 
weighted MP analyses of 18S + 28S rDNA 
and EF-la + Pol II sequences strongly and 
consistently corroborated the monophyly of 
Palpatores (Fig. 2), the result favored by anal¬ 
ysis of EF-la 4- Pol II alone, not the Dyspnoi 
4- Laniatores hypothesis that tended to be re¬ 
covered by the rDNA data analyzed alone (Gi- 
ribet et al. 1999). The strict consensus of 360 
trees (length, 6067; Cl, 0.31; RI, 0.48) derived 
from analysis of ribosomal and protein-encod¬ 
ing DNA included Opiliones, Phalangida, 
Laniatores, Palpatores, Dyspnoi and all rep¬ 
resented opilion superfamilies (except Travu¬ 
nioidea) but did not recover Eupnoi (Fig. 2A). 
The Dyspnoi 4- Laniatores hypothesis was re¬ 
covered in only 5% of bootstrap replicates and 
required 30 additional steps for recovery by 
parsimony analysis. Strict consensus of 78 
trees (length, 865; Cl, 0.54; RI, 0.75), derived 
from MP analysis of ribosomal DNA and EF- 
la + Pol II amino acids included the same 
major opilion clades listed above but also re¬ 
covered Eupnoi (Fig. 2B). The Dyspnoi + 
Laniatores hypothesis was recovered in only 
5% of bootstrap replicates and required 10 ad¬ 
ditional steps for recovery by parsimony. 

DISCUSSION 

Status and future of higher relationships 
in Opiliones.—Recent phylogenetic analyses 
clearly indicate that Opiliones consists of two 
clades, Cyphophthalmi and Phalangida (Giri- 
bet et al. 1999; Shultz 1998; Shultz & Regier 
unpubl. data). Results from analysis of elon¬ 
gation factor-la (EF-la) and RNA polymer¬ 
ase II (Pol II), both alone and combined with 
18S and 28S rDNA, strongly support the 
monophyly of Palpatores and Laniatores and 
are inconsistent with the recently proposed 
Dyspnoi 4- Laniatores hypothesis (Giribet et 
al. 1999, Giribet & Wheeler 1999). Molecular 
data examined thus far have consistently re¬ 
covered three palpatorean superfamilies (Cad¬ 
doidea, Phalangioidea, Troguloidea) as mono¬ 
phyletic groups (Giribet et al. 1999; present 
study), but monophyly of the remaining su¬ 
perfamily, Ischyropsalidoidea, is more prob¬ 
lematic. This superfamily was recovered by 
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Figure 2.—Strict consensus trees derived from combined, unweighted maximum parsimony analyses of 
18S rDNA, 28S rDNA, EF-la and Pol II. A. Strict consensus of 360 trees derived from analysis of rDNA 
and EF-la + Pol II nucleotides. R. Strict consensus of 78 trees derived from analysis of rDNA and EF- 

la + Pol II amino acids. Numbers above internal branches are bootstrap percentages. 


EF-la and Pol II in only a few analyses with 
low bootstrap support (Fig. 1) and in both 
combined analyses with moderate support 
(Fig. 2). The traditional grouping of Caddoi- 
dea and Phalangioidea into Eupnoi and Ischy- 
ropsalioidea and Troguloidea into Dyspnoi 
could be neither falsified nor confirmed by 
EF-la and Pol II, but they are both strongly 
supported in the combined analysis of 18S and 
28S rDNA conducted by Giribet et al. (1999). 


Rased on the present study, it appears that 
phylogenetic signal within EF-la and Pol II 
can resolve the deepest nodes within opilion 
phylogeny but become less useful at lower 
taxonomic levels. Specifically, relationships 
among and within superfamilies are largely 
unresolved by these two genes. This suggests 
that progress in opilion systematics will re¬ 
quire the development of new genes, espe¬ 
cially those that have evolved at a greater rate 
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than the two examined here as well as re¬ 
newed efforts to search for informative mor¬ 
phological characters. Expanded taxon sam¬ 
pling is also important, and particular 
emphasis should be placed on sampling rep¬ 
resentatives from basally divergent lineages of 
each well-defined opilion clade. These in¬ 
clude, but are not limited to, Crosbycus from 
Ischyropsalidoidea, neopilionids from Phal- 
angioidea, acropsopilionids from Caddoidea 
and the travunioid families from Laniatores. 

Relationships within Ischyropsalidoi¬ 
dea.—Ischyropsalidoidea is a morphological¬ 
ly diverse group encompassing three families 
(Xschyropsalididae, Sabaconidae, Ceratolas- 
matidae), but the monophyly of the superfam¬ 
ily is not well established. Shear (1986) united 
the ischyropsalidoids with four synapomor- 
phies, namely, metapeltidial cones, unseg¬ 
mented ovipositor, reduced palpal claw and 
presence of male cheliceral glands. However, 
metapeltidial cones are also present in Cad¬ 
doidea {Caddo) and the remaining characters 
occur in some or all members of Troguloidea 
(Shultz 1998). Further, our molecule-based re¬ 
sults do not support Shear’s (1986) system of 
ischyropsalidoid families (Fig. 3A), but nei¬ 
ther do they strongly conflict with it. Specif¬ 
ically, our analyses consistently recovered an 
Ischyropsalis + Ceratolasma grouping and 
frequently but weakly placed the ceratolas- 
matid Hesperonemastoma as the sister group 
to one or both sabaconids (Fig. 1). These re¬ 
sults suggest that Ceratolasmatidae and Sa¬ 
baconidae are not natural groups, although the 
evidence is weak. Still, para- or polyphyly of 
the two families is a hypothesis to be tested 
and is compatible with several other lines of 
evidence. First, few morphological characters 
support monophyly of the morphologically di¬ 
verse Ceratolasmatidae (= Acuclavella , Cer¬ 
atolasma , Crosbycus , Hesperonemastoma). 
Shear’s (1986) diagnosis of the family in¬ 
cludes no unique features and most characters 
have multiple states within the family (e.g., 
44 . . . pairs of tubercles on scutum . . . high and 
acute, or blunt and appressed ... or absent,” 
44 chelicerae with or without glands in males”; 
“palpi long, with ... or short, without. . . plu¬ 
mose setae”) (p. 13). Second, monophyly of 
Sabaconidae can also be questioned, although 
it is substantially more convincing than Cer¬ 
atolasmatidae. Potential sabaconid synapo- 
morphies include deep invagination in anterior 


midline of carapace, reduced sclerotization, 
and enlarged palpal tibia and tarsus (Shear 
1986). The first of these appears to be unique 
to the family, but the second seems to depend 
on a priori suppositions of character transfor¬ 
mation. The third character is undoubtedly de¬ 
rived and clearly expressed in Sabacon , but is 
less obvious in Taracus , especially when the 
palps are compared to those of Hesperone¬ 
mastoma (original observations). Thus, Saba¬ 
conidae is a probable but not unambiguously 
demonstrated family. Third, in his re-descrip- 
tion of Ceratolasma , Gruber (1978) noted two 
phenetic groupings of taxa within Ischyrop¬ 
salidoidea that are broadly congruent with our 
findings. Specifically, Ceratolasma and Ischy¬ 
ropsalis share a 4 ‘prominent sternum, large la¬ 
bium, palpi without plumose setae, but with 
numerous microtrichia, and also a complex 
midgut anatomy” (p. 109). He also noted that 
Hesperonemastoma is similar to sabaconids in 
having “less developed sterna, small labia, 
palpi with extensive development of plumose 
setae and reduction of the microtrichial cov¬ 
er” and simpler midgut anatomy. Admittedly, 
both lists are mosaics of primitive and derived 
traits, and more intensive morphological and 
molecular analyses are needed to make pro¬ 
gress in ischyropsalidoid systernatics. The 
molecular data are open to criticism in that 
relationships among the relevant taxa are un¬ 
stable and Hesperonemastoma appears to have 
undergone more rapid molecular evolution 
than other ischyropsalidoids, which may ac¬ 
count for the ambiguity. 

Relationships within Troguloidea*—In 
contrast to ischyropsalidoids, the troguloids 
are a well-defined, monophyletic group. Mar¬ 
tens (Martens 1980, 1986; Martens et al. 
1981; Martens & Suzuki 1966) and Shear & 
Gruber (1983) proposed several synapomo- 
phies, including a penis with two longitudinal 
muscles, unique unsegmented ovipositor, fu¬ 
sion of sternum and leg coxae, clavate palpal 
setae, and reduced palpal claws. Only the lat¬ 
ter is open to question, as reduced palpal 
claws are also present in ischyropsalidoids. 
The present study included one representative 
from each of the four troguloid families (Nip- 
ponopsalididae, Nemastomatidae, Dicranolas- 
matidae, Trogulidae) and our results strongly 
supported the monophyly of the superfamily 
under all character partitions and analytical 
methods (Fig. 1). However, relationships with- 
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in Troguloidea were ambiguous. Shear & 
Gruber (1983) regarded Dicranolasmatidae 
and Nemastomatidae as sister groups on the 
basis of one character (penis muscles with 
long tendons) but did not propose relation¬ 
ships between this clade, Trogulidae and Nip- 
ponopsalididae. Shultz (1998) proposed di- 
cranolasmatids and trogulids as sister groups 
based on heavy sclerotization and anteriorly 
projecting eye tubercle or “hood” in these 
taxa. Our molecular data do not strongly sup¬ 
port any phylogenetic arrangement among the 
troguloid families, although there was a ten¬ 
dency to recover the nemastomatid “Dendro- 
lasma” as the sister group to the remaining 
representatives. 
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